Photonic Hall effect in cold atomic clouds 
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Abstract: On the basis of exact numerical simulations and analytical calculations, we describe 
qualitatively and quantitatively the interference processes at the origin of the photonic Hall effect 
for resonant Rayleigh (point-dipole) scatterers in a magnetic field. For resonant incoming light, 
the induced giant magneto-optical effects result in relative Hall currents in the percent range, three 
orders of magnitude larger than with classical scatterers. This suggests that the observation of the 
photonic Hall effect in cold atomic vapors is within experimental reach. 
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Light propagation in homogeneous media in the pres- 
ence of a static magnetic field is a rich and vivid field of 
research where the symmetries dictated by nature lead 
to subtle magneto-optical phenomena [1 . About ten 
years ago, the question of magneto-optics in strongly 
scattering media was addressed and several effects bear- 
ing close analogies with electronic transport were theo- 
retically predicted and observed [2]. One striking exam- 
ple is given by the so-called photonic Hall effect (PHE) 
where light propagating in a scattering medium subject 
to a transverse magnetic field can be deflected in the 
direction perpendicular to both the incident beam and 
the magnetic field O |4] . Cold atomic gases provide an 
appealing testing ground for such effects in the multi- 
ple scattering regime. Indeed they constitute a perfect 
monodisperse sample of highly-resonant point scatterers 
which are very sensitive to external perturbations and 
where spurious phase-breaking mechanisms can be easily 
circumvented. Typically few Gauss are enough to induce 
strong magneto-optical effects like the Faraday rotation 
[5] in sharp contrast with classical materials where Teslas 
are needed to induce significant effects. If the impact of 
a magnetic field on coherent backscattering has already 
been studied [6l [7] , the question of the observation of the 
PHE in atomic vapors is still open. In this paper, we 
present analytical and numerical calculations identifying 
the physical origin of the PHE for point-dipole scatterers 
without internal degeneracies. Our results show that the 
effect should be observable in cold atomic gases. 

For a quantitative study of the PHE, one needs to ad- 
dress the question of directional asymmetries displayed 
by the configuration-averaged radiation pattern of an as- 
sembly of atoms located at random positions and illumi- 
nated by an incident monochromatic plane wave (wave 
vector k, polarization vector e ± k, angular frequency 
uj = ck = 27rc/A) while being subjected to an exter- 
nal static magnetic field with strength B pointing in 
the direction B. We consider here the simplest possible 



atomic internal structure serving our purposes, namely 
two-level atoms having a groundstate with angular mo- 
mentum J = connected by a narrow optical dipole 
transition to an excited state with angular momentum 
Je = 1. The energy separation between the atomic states 
is hwo = hc/Xo and the natural energy width of the ex- 
cited state is hT <C hwo. This is one of the best pos- 
sible natural realizations of resonant point scatterers |8] 
and it corresponds for example to the case of ^^Sr atoms 
(Ao = 461nm, r/27r = 32MHz, Lande factor of the ex- 
cited state ge = 1). When 5 = 0, the incident light is 
quasi-resonant (A ~ Aq) with this optical dipole transi- 
tion and we will denote hy 5 = {co — coq) the light de- 
tuning with respect to the atomic line {6 <C ujq). In 
the rest of the paper, we assume that the incident light 
intensity is low enough to neglect all nonlinear effects. 
When the magnetic field is applied, the internal degener- 
acy is lifted (Zeeman effect) and the excited level is split 
into 3 components separated by fiB where /i/27r = 1.4 
MHz/G is the Zeeman shift rate. As soon as the Zee- 
man shift becomes comparable to the resonance width, i.e 
(pB = ^jiB /T ^ 1, the scattering properties of each atom 
are strongly modified (this would occur at 5 ~ 11 G in 
the case of ^^Sr). 

The source of the field radiated by the atom is the 
oscillating electric dipole moment d exp(— icc;t) induced 
by the incident electric field Eexp(— zcjt). The radiated 
spectrum is elastic because there is no Zeeman effect in 
the groundstate. The situation would be more involved 
for atoms with a degenerate groundstate where frequency 
changes are possible, leading to an inelastic spectrum. In 
our situation, d = eoa(B)E and the radiation proper- 
ties are fully characterized by the polarizability tensor 
a(B) = ao ^(B) given by: 
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FIG. 1: (Color online) The photonic Hall geometry. A plane 
wave k = /cx is scattered by a cloud of atoms subjected to a 
static magnetic field B = Bz. The Hall current is measured 
either in the linear polarization channel y ^ x or in the 
y ^ z one. It is defined as AI = (/+ — /_), with I± the 
configuration- averaged differential cross- section along zby. 



where ao is the complex atomic polarizability at 5 = 0. 
The dyadic tensor T(B) embodies the effect of the mag- 
netic field on the photon polarization degrees of free- 
dom and gives rise to the usual magneto-optical effects. 
The C term is responsible for the normal extinction of 
the forward beam in a scattering medium made of such 
atoms (Lambert-Beer law). The rj term describes the 
magnetically-induced rotation of the atomic dipole mo- 
ment (Hanle effect) ^ and induces Faraday rotation and 
dichroism effects in the forward beam when k || B [1 . 
Finally, the ^ term is responsible for the Cotton-Mouton 
effect also observed in the forward beam when k ± B [T . 
These coefficients read: 
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and are real on resonance ((^ = 0). From the polarizabil- 
ity we then get the single-atom differential cross-section: 
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where ctq = |aop/c^/67r is the total scattering cross- 
section at zero magnetic field. As immediately seen, the 
single atom differential cross-section ([3| only depends on 
the incoming and outgoing polarizations and is thus com- 
pletely insensitive to the change k' — k^ For an iso- 
lated atom there is thus no possible asymmetry when 
reversing the direction of observation [9l [10] . As a con- 
sequence, the single scattering signal originating from an 
assembly of such atoms cannot display any directional 
asymmetry and the PHE, if any, must come from a mul- 
tiple scattering effect. 

Before considering the general case, we first analyze the 
radiation properties of two isolated atoms located at po- 
sitions ri and r2. This is the simplest possible situation 
where multiple scattering plays a role and it will allow us 
to extract physical insights about possible mechanisms 
at work in the PHE. In the Hall geometry (see Fig. [T]), 
one measures the differential cross sections I± along the 
up and down directions ±y perpendicular both to the 



incident light direction k = /cx and to the magnetic field 
B = 5z. The total radiation field is the sum of the fields 
radiated by the atomic dipoles induced by the incoming 
and scattered fields at their respective positions. The ex- 
act solution thus involves the inversion of a linear system 
of 2 coupled vectorial equations where the polarizabil- 
ity tensor plays a key role. For a fixed relative distance 
r = |ri — r2| between the atoms, we compute exactly the 
differential cross sections I± and average them over all 
possible relative orientations of the atoms to get (/±). 
We then extract the Hall current A/ = — (/-), the 
mean intensity / = ((/+) + {!-))/ 2 and the relative Hall 
current P = AI/I. There are 4 possible linear polar- 
ization channels e ^ e' for the data analysis. However 
we first note (and we have numerically checked) that the 
Hall currents must be the same in the linear channels 
y ^ z and z ^ x since they are related by time-reversal 
symmetry. Second, as we also checked numerically, the 
Hall current must vanish in the z || z channel since the 
polarizations being along the magnetic field, they are in- 
sensitive to the Hanle effect. We are thus left with the 
two lin -L lin channels y ^ x and y ^ z. Fig. [2] summa- 
rizes our numerical results. We see that the Hall current 
vanishes as kr 0. Indeed for very small distances, the 
two radiating dipoles are always in phase and thus add 
up constructively meaning that the two atoms behave 
like a single scatterer with a dipole moment twice larger, 
a situation for which we already know there is no pos- 
sible directional asymmetry. We also see that for low (/) 
values, the relative Hall current decreases in the y ^ z 
channel, whereas, in the y ^ x channel, it is comparable 
to the one for high values. The reason is that, in the 
y ^ X channel, the single scattering due to the Hanle 
effect increases with the magnetic field, such that the 
background intensity in this channel gets more and more 
contaminated (and even dominated at large distances) 
by the single scattering signal. In the y ^ z channel, it 
is always filtered out. 

To get some insights about the physical processes at 
work, we consider the case of a "dilute" medium kr ^ 1 
and we expand the field radiated by the two atoms in 
powers of {kr)~^. Skipping tedious details, at leading or- 
der the scattered amplitude is obtained from the sum of 
the two diagrams shown in Fig. |3^, the respective ampli- 
tude being denoted by u and v. Each differential cross- 
section I± contains thus interference terms (i.e. uv-\-vu) 
and background terms (i.e. Ii^p + I'up), but since the later 
are independent of the scattering direction, they cancel 
out in the difference. The Hall current is thus solely given 
by a difference of interference terms, which are precisely 
the crossed terms at the heart of the coherent backscat- 
tering effect [TIJ. More precisely, the Hall current reads 

A/= ((5/(r,B){cos[(k + kO -r] -cos[(k-kO -r]}), 

(4) 
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FIG. 2: (Color online) Two- atom case. Relative Hall current 
(3 observed at resonance ^ = in the linear polarization chan- 
nels y ^ X (a) and y ^ z (b) as a function of the relative 
distance kr for different values of = 2/iB/r. 



average over the relative orientation of the two atoms 
and Ar is the projector onto the plane perpendicular to 
r. After elementary manipulations, equation Q can be 
conveniently rewritten as follows: 

A/ = ( ((5/(r,B) -(5/(r,-B)) cos[(k + • r]) . (5) 

The Hall current can then be understood as a difference 
between the two configuration-averaged interference ef- 
fects generated in the same direction k', but for opposite 
directions of the magnetic field. The imbalance results 
from the interplay between the dipole rotation induced by 
the Hanle effect and the transverse projector Ap. More- 
over, since the cosine term depicts very fast oscillations, a 
stationary phase approximation shows that the main con- 
tribution in Eq. ([5| comes from the configurations where 
the two atoms are aligned along the k + k' direction. 

Performing the angular average, one gets the Hall cur- 
rent in the y ^ z channel: 
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the \/2kr term simply corresponding to the value of (k + 
k') • r when these two vectors are parallel to each other. 

In the y ^ X channel, a closer inspection shows that, 
at same order, there is an additional contribution due to 
diagrams accounting for the interference between recur- 
rent and single scattering processes (see Fig. [s})). The 
final expression for the Hall current in this channel is 
quite tedious but simplifies at resonance (5 = 0: 
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As shown in Fig.|4]for (5 = and (/)^ = 1, the agreement 
between the numerical computations and these asymp- 
totic results proves excellent as soon as /cr > 10. In 



FIG. 3: (Color online) Two-atom case. For large distances, 
the PHE in the y ^ z channel results from the interference 
between the two diagrams (a). In the y ^ x channel an 
additional contribution comes from the interference between 
the two diagrams accounting for recurrent scattering (b) (and 
also between the ones obtained by exchanging the two atoms) . 
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FIG. 4: (Color online). Normalized Hall current in the 
two-atom case. Comparison between the numerical calcula- 
tions (plain curves) and the analytical predictions Eqs. (|6]),([7| 
(dashed curves) as a function of kr for ^ = and 0^ = 1. (a) 
y ^ z channel, (b) y ^ x channel. The excellent agreement 
emphasizes that the underlying physics of the PHE is fully 
catched by the diagrams depicted in Fig. |3] 



particular it is quite remarkable that recurrent scatter- 
ing is essential to reproduce the additional oscillations 
observed in the y ^ x channel and due to the cos2/cr 
term in Eq. ([7|). We have also numerically checked that 
the agreement is still excellent in both channels when 
(5 7^ 0. This gives clear evidence that the physical mecha- 
nism at the heart of the PHE for atomic scatterers is the 
interference between the scattering processes depicted by 
the diagrams in Fig. |3] In addition, the fact that, in the 
small (\)B limit, the Hall current is proportional to (/)^ (i.e. 
to the atomic dipole rotation) actually reflects that only 
one active scatterer would be enough to generate a Hall 
current. Finally, one must also mention that, contrary 
to the case studied in Ref. [3 , both our analytical and 
numerical results show that there is still a PHE when the 
antisymmetric part of the self-energy (S oc a(B)) has a 
vanishing real part (i.e. at resonance (5 = 0), the exact 
dependence on (5 being a smooth bell-shaped curve. 
Assuming an isotropic propagation in the medium 
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TABLE I: The left side displays the relative Hall currents 
P at 6 = and = 1, in the Hall geometry, for 500 atoms 
uniformly distributed inside a sphere with optical density p ~ 
0.65. The number of disorder configurations is 6 X 10^ The 
expected accuracy of the order of 10~^ is confirmed when 
comparing to the relative currents p (see right side) in the 
k II B configuration, when no PHE should show up. 
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(with a mean free path one can further average the 
preceding expressions over all possible positions of the 
two scatterers. As seen from Eqs. ([6| and ([7|, the in- 
terference effect leads to fast decreasing oscillations at 
the wavelength scale. Since A <C ^, the main contribu- 
tion to the Hall effect will arise for scatterers separated 
by r ~ A. The probability of finding two such scatter- 
ers scales with the optical density p = nA^, n being the 
scatterer density; we thus expect the actual relative Hall 
current /3 to be smaller by a factor {ki)~^ than the av- 
erage background [12]. Obviously, the exact result will 
depend on the geometry of the medium and on the op- 
tical thickness. Nevertheless, the PHE in a assembly of 
resonant point scatterers, even if rather small, should be 
measurable. 

To bring a numerical confirmation of our findings, we 
have considered N = 500 atoms uniformly distributed 
inside a sphere at an optical density p ^ 0.65 and illu- 
minated by a plane wave set on resonance. This leads 
to M = 47rV3p ^ 20 at 5 = 0. Such a value is difficult 
to achieve in a real experiment but has already been ob- 
tained [13] . The optical thickness along a diameter of the 
sphere isb ^ 3.5. The magnetic field value has been set at 
= 1- To obtain the total radiated field, we have solved 
the corresponding system of 3A^ linear equations [M] and 
we have computed the various quantities of interest av- 
eraged over 6 X 10^ configurations, such that we expect 
an accuracy of the order of 10~^. To stress the existence 
of the PHE, we have also computed (3 in the geometry 
k II B where no PHE should be observed; as depicted 
by the right side of table [l| the corresponding values are 
effectively at most of the order of few 10~^. Up to this 
accuracy, the numerical results in the geometry k ± B, 
are in a qualitative agreement with the two- atom case: 
there is no Hall current in the z ^ z channel. It is about 
the same order of magnitude in the two conjugate chan- 
nels y ^ z and z ^ x (/? 5.5%). Finally, it is larger 
in these channels than in the y ^ x channel (|/?| ~ 1%). 
To enforce the validity of our mesoscopic description, we 
have computed /3 for different values of ki^ while keeping 
fixed the optical thickness b. The results are displayed in 
table [H] and, as expected, the product pk£ is, within the 



TABLE II: Relative Hall current f3 in the z ^ x and y ^ z 
channels at a fixed optical thickness 6 = 1.3. Within the nu- 
merical accuracy, the results indicate that p scales like l/{ki). 
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statistical errors, independent of ki. Finally, one must 
note that the values found here are in the percent range 
and are much larger (by three orders of magnitude) than 
observed with classical scatterers [4 although the mag- 
netic field is smaller by two or three orders of magnitude. 
This arises because the Zeeman effect in highly-resonant 
atomic scatterers induces a "giant" dipole rotation which 
amplifies the PHE. 

In summary, on the basis of numerical and analytical 
calculation, we have qualitatively and quantitatively ex- 
plained the underlying interference effect at the origin of 
the photonic Hall effect for resonant point-dipole scatter- 
ers. Our results show that the effect, albeit small, should 
be observable in cold atomic vapors. Possible extensions 
of this work would consist first in developing the dia- 
grammatic analysis to an arbitrary number of scattering 
events and, second, in accounting for internal degenera- 
cies in the atomic groundstate. This would give a quan- 
titative and comprehensive description of the photonic 
Hall effect in cold atomic clouds. 

The Authors thank G.L.J. A. Rikken and 
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